Dissipative effect of thermal radiation loss in high-temperature dense plasmas 
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A dynamical model based on the two-fluid dynamical equations with energy generation and loss 
is obtained and used to investigate the self-generated magnetic fields in high-temperature dense 
plasmas such as the solar core. The self-generation of magnetic fields might be looked at as a self- 
organization-type behavior of stochastic thermal radiation fields, as expected for an open dissipative 
system according to Prigogine's theory of dissipative structures. 
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Thermal motion of electrons generates thermal radia- 
tion mainly by means of thermal bremsstrahhmg emis- 
sion in high-temperature plasmas. The higher the plasma 
temperature the more intense thermal radiation of the 
plasma. Thermal radiation causes substantial energy loss 
in high-temperature plasmas such as thermal nuclear fu- 
sion plasmas; for instance, the well-known Lawson condi- 
tion is derived from the condition that the fusion energy 
compemsates the thermal bremsstrahlung emission loss. 
Here we are about to report some dissipative effect of 
thermal radiation loss in high-temperature dense plas- 
mas such as thermal nuclear fusion plasmas. 

Generally, the elastic collision frequency u e i between 
electrons and ion is high in high-temperature dense plas- 
mas, for example, in the core of the Sun T e ~ 10 7 K, 
n e ~ 10 32 to~ 3 , we can eastimate v e i ~ 10 16 Hz, the 
plasma frequency v pe ~ 10 17 Hz. Such a frequent colli- 
sion would destroy any collective plasma motion or co- 
herent structure if the Sun were isolated adiabatically. In 
other word, the Sun would be in thermodynamic equilib- 
rium if it was an isolated system. However, the Sun is 
in fact an open system: it not only gains thermal en- 
ergy from the inelastic collisions between ions such as 
the proton-proton chain (or pp chain), but also loses en- 
ergy by means of thermal radiation, for example, the 
loss rate of thermal radiation or the photon luminosity 
L Q ~ 4 • 10 33 erg/s (Bahcall, 1989). Therefore, the Sun 
may be far from equilibrium. The theory of dissipative 
structures (Prigogine, 1973; 1977) has shown that co- 
herent behaviors may occur in a dissipative system. In 
fact, we have witnessed and are witnessing the dissipa- 
tive structures in the convectve zone of the Sun, the con- 
vective cells. These are the coherent structures at the 
outer part of the Sun caused by thermal radiation dissi- 
pation which maintains the necessary temperature gra- 
dient. What are the dissipative structures in the source 
region? 

Roughly, nuclear fusion reactions take place within 
0.3i?© (Bahcall, 1989), where nuclear reactions (inelastic 
collisions between ions) first enhance stochastic motion 
of ions, the ions then transfer energy and momentum to 
electrons via elastic collision. The electrons lose energy 
and momentum through emitting photons via inelastic 



collision with ions and the thermal radiation escapes from 
the source region via radiation transfer. On average, the 
direction of the energy and momentum flux in the nuclear 
fusion region is therefore: 



ions — > electrons — > photons 



(1) 



This implies that electrons do not reach thermodynamic 
equilibrium with ions (Tj > T e ), nor do photons with 
electrons in the fusion plasma and that it is the very 
ions that replenish the radiation loss of both energy and 
momentum of electrons. 

This noncquilibrium effect also shows itself in the two- 
fluid dynamical equations (Rose and Clark, 1961; Ma et 
al, 1988; Rybicki and Lightman, 1979): 
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where ds is a differential element of length along the ray, 
I v is the energy crossing unit area in unit time, unit 
solid angle and unit frequency range of radiation fields, 
called as the specific intensity or brightness, a„ and j v 
are the absorption coefficient and emission coefficient re- 
spectively. P a is called the kinetic stress tensor, the 
thermal flux vector, F Q is the total force exerted on the 
fluid element of a species, Of five collision-related terms, 
8n a /St represents created or destroyed particle density of 
a species in the unit time, R Q is the momentum gain or 
loss rate of the fluid element of a particles by their elas- 
tic collision with the unlike particles (including photons) , 
the same for (5R a but due to the inelastic collision. Sim- 
ilarly, Q a and SQ a are thermal energy gain or loss rate 
due to the elastic collision and the inelastic collision re- 
spectively. 



1 



For the sake of clarity, we assume that the quasi- 
stationary state has been reached with dT a /dt w 0, as 
believed in the solar core. If we further neglect viscosity, 
thermal conduction, thermal convection and restrict to a 
very short spacescale (< 1cm), the energy conservation 
equation reduces to 



= Q a + 5Q a 



(6) 



Eq.(g) tells us that radiation energy is continuously 
transfcrcd from the nuclear fusion region to the outside 
due to the radiation loss, which implies that electrons 
lose energy through emitting photon via inelastic collision 
with ions, SQ e < 0, Eq.(||) thus leads to Q e = —5Q e > 0. 
This means that electrons should gain energy on average 
from ions via elastic collision with ions. In the quasi- 
stationary case, Qi = J2 a ^ e Q a ~ ~Qe < 0, which im- 
plies that ions lose energy on average via elsatic collisions 
with electrons. The energy loss of ions Qi < will be re- 
plenished by the energy generation 5Qi > due to the 
inelastic collision (i.e., nuclear reaction) between ions or 
ions and electrons, Qi = —8Qi, as implied by the energy 
conservation equations for ions. The energy transfer ac- 
companies similar momentum transfer between electrons 
and ions: R e = — Ri = ^2 a ^ e Ra > 0, which means that 
ions lose momentum and electrons gain momentum on 
average in the elastic collision between ions and electrons 
in the quasi-stationary case. The momentum generation 
<5Ri = J2 a ^e 3R-a > of ions due to the nuclear fusion 
collision replenishes the momentum loss R^ < of ions. 

The following task is to estimate R^ + JRj and R e + 
<5R e . In so doing we have to take into account the radi- 
ation fields. The transfer equation (Rybicki and Light- 
man, 1979) takes a particularly simple form if, instead of 
s, we use the optical depth r v defined by dr u = a v ds, 



dlv 
dT v 



-Iu +j v /a v 



(J) 



In plane-parallel media, a standard optical depth is some- 
times used to measure distance normal to the surface, so 
that r„ = t v (z). A medium is said to be optically thick or 
opaque when t„, integrated along a typical path through 
the medium, satisfies t v > 1. An optically thick medium 
is one in which the average photon of frequency v can- 
not traverse the entire medium without being absorbed. 
Solar core plasma is optically thick, and the absorption 
is mainly due to the (massive) ions (Bacall, 1989). If a 
medium absorbs radiation, then the radiation exerts a 
force on the medium, because radiation carries momen- 
tum. The radiation pressure F ra( j can be defined by the 
specific intensity /„ in remembering that a photon has 
momentum hv/c, 
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, ,../ = - / a v I v ndQ,dv 



(8) 



z-direction, for instance. Since the absorption is mainly 
due to the ions, we know 



Rz + <5Ri = F r 



R e + SR e « 



(9) 
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As a result, the momentum transport equation reduces 
to 



mini—— = - Vpi + n l F l + F rad 
dt 



dV e 

1 dt 



= -Vp e + n e F e 



(11) 



(12) 



where to,, rii and Vj are the effective mass, density of 
ions and the effective velocity of an ion-fluid element. 
This is in agreement with the nonequilibrium processes 
(0) . The most outstanding characteristic of the reduced 
equation is that the collision terms have been canceled 
out. In contrast, the cancellation cannot occur in the 
equilibrium case because there is no a well-determined 
direction of energy and momentum flux. Consequently, 
we find out that the nonequilibrium effect cancels the 
collision effect. 

In the two-fluid approximation V Q = V Q (r, t) is the 
velocity of an a fluid element, which is the first-order 
moment of the distribution function at the velocity space 
of the a particle. Under such an approximation, the ac- 
celeration behavior of single particles has been ruled out. 
Therefore, the Maxwell equations (Rose and Clark, 1961; 
Ma et al, 1988) 



M(E, B; p, J) = 



(13) 



with induced charge density p = ^ a q a n a and induced 
current density J = ^ a n„(j Q V a determine only the in- 
duced fields in the plasma by motion of fluid elements. 
The radiation fields E rac i and B ra( j are generated by the 
acceleration motion of single particles and thus can be de- 
termined only by the retarded potentials of single moving 
charges, or Lienard-Wiechart potentials. On averaging 
single particle fields E S i„ and B S i„ in the velocity space 
with the distribution function of electrons as the weight 
factor, we can obtain the radiation fields (Rybicki and 
Lightman, 1979). At this stage, we are able to write 
down explicitly 

[(E + E rad ) + V a x (B + B rad )] + g Q (14) 

where g Q is gravity, a = i, e. The quasi-stationary con- 
dition requires that the thermal and radiation pressure 
balance the gravity on the macroscopic level (3> 1cm), 



■ft = o 



(15) 



where n is a unit vector along the direction of the ray, 



Vp e + g e = 



(16) 
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As a result, Eqs.fll l|) and (tt2h reduce to 

m a n a ^ = -J a T a Vn a + n a q a (E tot + V Q x B tot ) 
at 

(17) 

where h a (a = i, e) is the microscopic density fluctuation 
of the a fluid. E tot = E + E rad , the same for B tot . We 
have assumed thermal pressure p a = 7 Q n tt r Q with the 
polytropic index j a . The mass consernvation equation 
can be rewritten down as follows by neglecting particle 
creation or destruction: 



~dt 



V • (n a V a ) = 



(18) 



So far we have shown that the plasma in the nuclear fu- 
sion region can approximately be described by Eqs.(p_3|), 
( |l7| ) and ( ]l8| ) when the energy and momentum flux (Pis 
well-defined by the thermal radiation loss. These equa- 
tions possess two characteristics: (a) The elastic collision 
terms have been canceled out by the inelastic collision 
terms; (b) thermal radiation fields only appear in the 
momentum transport equations. The first confirms the 
fact that non-equilibrium may be a source of order, the 
second allows us to study the self-organization phenom- 
ena (Nicolis and Prigogine, 1977) of the random thermal 
radiation fields. Self-organization will lead to order, or 
dissipative structures. 

Self-organization needs nonlinearity. The reduced two- 
fluid dynamical equations contain the needed nonlin- 
earity. These equations can be simplified substantially 
through distinguishing the ion-timescale Tj ~ w" 1 and 
the electron-timescale r e ~ us~} with the self-generated 
effect of magnetic fields included in (Kono et al, 1981; Li, 
1993a,b): 
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where E is the slowly varying complex amplitude of 
the induced high-frequency (uq) electric field E: E = 
i{E exp(— uiot) + c.c.}. The plasma has been assumed to 
be plane-parallel to the density gradient of scalelength 
L along z-axis. Therefore, all variables depend on only 
z. Then, k = (0,0, k). E* = [E x ,E y ] is the transverse 
field component, while E z is the longitudinal field com- 
ponent. The validity conditions of these equations are 
k/koe <C 1, Wk/k De <C 1 and \Sn\/n e < 5v/v pe <C 1. 
Another physical constraint is Sn/n e < 8v/v pe . All vari- 
ables have been rescaled as follows (Morale and Lee, 



1977): B = B S /B ML , E = E/E ML , n = {Sn/n e )/N M L, 
t = t/TML and z = z/Zml, in which the tilded stand for 
the unrescaled. B ML = (3ia 2 W, E ML — 2(3iaE , N ML = 
(3ia 2 W, Tml — 2a//3iujp e and Zml = cT X ^l\[Wi are 
the new units, in which a = (Luj pe /\/3VTe) 2 ^ 3 , Pi — 
l+3(Ti/T e ), W = E$/{n e KT e /e ). The other parameters 

ft 3 , p = a 3 W/(l 



are: a = (c/V Te ) 2 /3, (3 



Bl) in that 



B a = B s /(m e u) pe /e), 7 = fa 3 (m e /TO;), 7' = 7/(1 -Bl) 2 , 



S x = (l/2/?f)cos0sin0, 



Sy — 



(l/2/3f)sin^sin0 and 



S z = (l/2ft 2 ) cos 6 in that both 9 e [0, n] and <p G [0, 2ir] 
are random due to random orientation of the thermal 
radiation field. Eq sw [p(v pe , T)5v / e^] 1 / 2 is the thermal 
radiation field components at and near the local plasma 
frequency v pe , where p(v,T) is spectral energy density 
of thermal radiation. The symbol "*" means convolu- 
tion product, i>i is the Landau damping of low- frequency 
longitudinal waves. 

The so-called self-generated magnetic field B s of the 
plasma is the slowly varying component of the induced 
magnetic field (Kono et al, 1981): 



B. 

B r 



Vr eN2 iE x E * 



= ( — ) 



(22) 



where B c = m e oj pe /e is the critical magnetic field at 
which electron gyrofrequency f2 e = eB c /m e equals to lo- 
cal plasma (angular) frequency uj pe , E c — (n e nT / to) 1 / 2 . 
The magnetic-field generation was shown to be due to a 
solenoidal current j s — — ieWp e /(167rm e Wg)V x (E x E*) 
[Eq.(l.l) of Kono et al, 1981] by many authors (see Kono 
et al, 1981 and references cited therein). Through taking 
the curl of the solw-timescale Ampere law and neglecting 
the displacement current, one obtains [Eq.(4.7) of Kono 
et al, 1981] 



47T 

VxVxB s = — V x j s 

c 



(23) 



This equation reduces to Eq.(|22|) when V ■ B s = is sat- 
isfied, as one of the Maxwell equations requires. Eq.J^) 
shows that if there is no self-organization of the stochas- 
tic thermal radiation field E rac i, which has been incorpo- 
rated in parameters S and p, the self-generated magnetic 
field B s must be very weak, vice versa. Therefore, we 
can investigate the self-organization effect of the ther- 
mal radiation field through monitoring B s . In order to 
make sure that we have observed the self-organization 
phenomenon, we have to use the statistical average Bq(z) 
of B s (z',t) over time, space and random initial condi- 
tions: 



B (z) = 



1 



Bi — B rl 



-B r 



(24) 



where B t = L^/™"* [B 2 (z') + B 2 y (z')]dz'Y' 2 and 
n is the number of solutions starting from random ini- 
tial conditions. B means averaging B in the range 
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B max /B c € [0.8,0.9], where B max = Max{[B 2 x (z) + 
By(z)] 1 / 2 } < 0.9B C is one of our controlling parameters 
(the other is the maximum density fluctuation 8n/n e < 
Sis /is). L s i m is the simulation cell size. 
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FIG. 1. Variation of the self-generated magnetic field with 
the solar radius coordinate z/Rq, where Rq is solar radius. 
Curves I and 2, representing solar plasma temperature and 
electron number density respectively, are the input parame- 
ters. Symbols "+" stand for the corresponding self-generated 
magnetic fields. 
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A standard second-order explict quasispectral method 
(Li and Li, 1993; Li, 1996) has been used to numerically 
solve the model equations in order to calculate B x and 
By (note: B z = 0) and thus B s . We have used table 
XVI of Bahcall & Pinsonneault 's standard solar model 
(1992) to calculate our input parameters: the electron 
and ion temperatures T e as Tj = T, electron number 
density n e = ^(l+X)p/mu, and plasma scalelength L = 
n e (dn e /dz)^ 1 , where X = X^H) + X( 3 He) represents 
the concentration, by mass, of hydrogen. Although the 
simulation cells are small, L S i m ~ 10 -3 m, it is much 
larger than the Debye length \o e = Vre/wpe ~ 10 1 
m. Figure [I] shows Bq(z) obtained by using n = 10 (the 
convergence trend has been found by using n = 1,4,8 
and 10), which shows at least 0.4 gigagauss magnetic 
field may be generated at the center of the Sun. Such a 
strong field shows that the self-organization behavior of 
the stochastic thermal radiation fields does occur. The 
collision time scale is about 10~ 16 s, while the growth 
time scale of the self-generated magnetic field is about 
10~ 12 s. The fact that the coherent time scale is much 
longer than the collision time scale is also an indicator of 
self-organization (Nicolis and Prigogine, 1977). 
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